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Photonic crystals and metamaterials represent two seemingly different classes of artificial elec¬ 
tromagnetic media but often they are composed of similar structural elements arranged in periodic 
lattices. The important question is how to distinguish these two types of periodic photonic structures 
when their parameters, such as dielectric permittivity and lattice spacing, vary continuously. Here, 
we discuss transitions between photonic crystals and all-dielectric metamaterials and introduce the 
concept of a phase diagram and an order parameter for such structured materials, based on the 
physics of Mie and Bragg resonances. We show that a periodic photonic structure transforms into 
a metamaterial when the Mie gap opens up below the lowest Bragg bandgap where the homoge¬ 
nization approach can be justified and the effective permeability becomes negative. Our theoretical 
approach is confirmed by detailed microwave experiments for a metacrystal composed of a square 
lattice of glass tubes filled with heated water. This analysis yields deep insight into the properties 
of periodic photonic structures, and it also provides a useful tool for designing different classes of 
electromagnetic materials in a broad range of parameters. 


The concepts of phase transitions and phase dia¬ 
grams are among the most fundamental concepts in 
physics. Here we reveal and analyse a rather unusual 
type of phase transitions between two classes of artificial 
electromagnetic structures: photonic crystals (PhCs)i2 
and metamaterials (MMs^. Rapidly growing atten¬ 
tion to the study of PhCs and MMs is due to their 
unique electromagnetic properties and the exciting new 
physics these structures may offer for novel applications 
in photonics^ - —. The resonant properties of PhCs are 
driven by the familiar Bragg scattering from a periodic 
structure when the lattice spacing becomes comparable 
with the wavelength of light. MMs are regarded as ef¬ 
fective media with averaged parameters which in turn 
can be engineered by changing the properties of their 
subwavelength elements modifying the propagation of 
electromagnetic waves with the wavelengths much larger 
than the lattice spacing. 

Many properties of conventional MMs are usually 
defined by the optically-induced magnetic response of 
metallic wires or split-ring resonators, which has been 
demonstrated experimentally in many studies for both 
microwave and near-infrared spectral region s 10 ! 11 . How¬ 
ever, the scaling of such MMs to higher frequencies en¬ 
counters many problems, including strong absorption 
of electromagnetic radiation in metals. Recently sug¬ 
gested electromagnetic periodic structures composed of 
high-permittivity dielectric elements^ - — can avoid these 
problems, while still delivering an effective magnetic re¬ 
sponse. The main physics underlying the properties 
of such all-dielectric MMs is based on the resonant 
modes of the Mie scattering^ supported by dielectric 
particle s 16 ! 17 . Mie resonances provide a mechanism for 
the generation of an optically-induced resonant magnetic 
response based on the displacement current, and they 


offer a simpler and more versatile route towards the fab¬ 
rication of isotropic MMs and metasurfaces operating at 
optical frequencies. 

PhCs and all-dielectric MMs can be composed of the 
similar structural elements arranged in the similar ge¬ 
ometries. Such structures offer a unique opportunity to 
study a transition from PhCs to MMs when their param¬ 
eters, such as dielectric permittivity and lattice period, 
vary continuously. The transition we describe here can be 
mapped onto a phase diagram; however it is quite differ¬ 
ent from phase transitions well known in other fields such 
as thermodynamics. The thermodynamical phase tran¬ 
sitions appear as transformations of a system from one 
state or “phase” to another with changing external vari¬ 
ables such as pressure, temperature, chemical potential, 
etc. Here we reveal and analyze a rather unusual type 
of phase transitions between two regimes of light prop¬ 
agation, i.e. the regime where the propagation of light 
is mainly influenced by the interference between multiple 
Bragg scattered (PhC “phase”), and the regime where 
the propagation of light is mainly determined by the 
properties of each single element of the materials (MM 
“phase”). 

First, we identify two different scenarios for the transi¬ 
tion of a two-dimensional square lattice of dielectric rods 
from PhC to MM (see Fig. 1). The first scenario that has 
never been discussed previously occurs when the lattice 
spacing a decreases in comparison with the fixed radius 
of rods r and dielectric permittivity e. This transforma¬ 
tion leads to an increase of the filling ratio r/a, and it 
makes possible the homogenization^ of the periodic di¬ 
electric structure with the negative effective permeability 
(pL < 0) at higher values of £—— . 

The second scenario occurs when the dielectric per¬ 
mittivity of rods 5 varies in a fixed square lattice. This 
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FIG. 1: Transition from photonic crystals to metama¬ 
terials. Schematic illustrates two possible transformations 
of the PhC b to MM a,c through the competition between 
the Bragg and Mie scattering wavelengths. To satisfy the 
classical homogenization condition in the low-frequency long- 
wavelength limit (A Bragg < A Mie), one may decrease the lat¬ 
tice spacing and increase the filling ratio r/a (r/a- scenario) 
or, alternatively, increase the rod permittivity e (e-scenario). 


transformation leads to an increase of the Mie scatter¬ 
ing wavelength, and it can also provide the conditions 
for the homogenization approach with negative y. While 
photonic structures with varying e were studied earlier 
theoreticall y 22 ! 23 , the transition between PhCs and MMs 
has never been analyzed previously to our best knowl¬ 
edge, especially in terms of the phase diagram, as dis¬ 
cussed below. 


RESULTS 

The interplay between Bragg and Mie resonances 

To analyze the PhC - MM phase transition, we in¬ 
vestigate theoretically electromagnetic properties of the 
dielectric periodic medium depending on the dielectric 
contrast, filling ratio and incident wave polarization. As 
an illustrative example, we consider here 2D square lat¬ 
tices of parallel dielectric circular rods infinitely long in 
the z direction with the radius r and real frequency in¬ 
dependent permittivity 5. The rods are embedded in air 
(e a i r = 1) and the arrays can be characterized by the fill¬ 
ing ratio r/a. We have derived three key sets of spectro¬ 
scopic data (see Methods), namely: (i) spectra of the Mie 
scattering by an isolated rod; (ii) photonic band structure 
of an infinite 2D square lattice composed of rods; and (iii) 
transmittance of a 2D square lattice of 10 layers in length. 
All data sets have been calculated for a wide range of the 
rod permittivity 1^6^100 with the step of = 1 (e- 
scenario) and in the range of filling ratio 0^r/a^0.5 with 
the step of 0.01 (r/a- scenario). Note that at r/a = 0.5 all 
rods touch each other and at r/a > 0.5 the rods penetrate 
into each other and the structure appears as inverted one 


of air holes in a dielectric matrix. The calculated data 
(see Methods) allow us to analyze precisely the evolution 
of the low-frequency region of the photonic band struc¬ 
ture formed by a complicated mixture of the dispersion 
curves originating from Mie and Bragg resonances and 
propagating modes. Since our goal is to discuss the PhC 
- MM transition, we concentrate on the TE polarization 
where the lowest TEoi magnetic dipole resonance of the 
rods results in a negative magnetic permeability in the 
Mie gap spectral region^. 

The physics of the PhC - MM transition is displayed 
by the gap map (Fig. 2) that relates the Bragg and 
Mie gaps positions to the rod permittivity 5 for given 
r/a = 0.25. These illustrate essential features such as 
strong decreasing of Mie resonance frequencies with 5 
increase followed by intersection of the Mie and Bragg 
bands. The Mie resonances of individual rods form both 
non-dispersive flat bands in the photonic structure and 
Mie photonic gaps in the energy spectrum (Figs. 2,6). 
For low-contrast PhCs, all Mie resonances are located is 
a frequency range higher than the lower Bragg gap. With 
increasing rod permittivity e, the Mie resonances of sin¬ 
gle rods and correspondingly Mie bands of the structure 
demonstrate a shift from the higher to low frequencies, 
crossing the strongly dispersive Bragg bands. The cre¬ 
ation of the lowest Mie gap driving the artificial mag¬ 
netism is clearly seen both in the band structure and in 
the transmission spectra (Figs. 2,6). 


A metacrystal: two crossover scenarios 

Armed with the results of our calculations, we verify 
our concept experimentally by engineering a metacrystal 
composed of plastic circular tubes filled with water and 
forming a 2D square lattice with variable lattice constant 
a (Fig. 3). The structure of the metacrystal enables us 
to perform two types of experiments according to the r/a 
- and ^-scenarios. The r/a -scenario was realized using a 
special construction of the metacrystal that makes it pos¬ 
sible to change the distance between tubes synchronously 
in x and y directions keeping the 2D square symmetry 
of the lattice (Fig. 3). For the ^-scenario, we employ 
the advantages of strong temperature dependence of the 
dielectric permittivity of water that is characterized by 
5 = 80 at 20°C and £ = 50 at 90° C in the microwave 
frequencies range from 1 GHz to 6 GHz2i. Based on the 
calculated results and the parameters of the metacrys¬ 
tal, we chose intervals of most interest for each scenario, 
namely 60^5^80 at 0.7^r/a^0.8 for the ^-scenario and 
0.65^r/a^l.9 at e = 60 for the r/a- scenario. 

The transmission spectra of the metacrystal installed 
in an anechoic chamber were measured in the 1 GHz to 
3 GHz range. A rectangular horn antenna (TRIM 0.75 
GHz to 18 GHz; DR) connected to a transmitting port 
of the vector network analyzer Agilent E8362C was used 
to approximate plane-wave excitation. A similar horn 
antenna was employed as a receiver. For the electromag- 
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FIG. 2: Bandgap diagram for the 2D square lattice of circular rods, a The bandgap diagram obtained from the 
photonic band structures calculated for 1^5^60 with steps of Ae = 1. The rods are embedded in air, r/a — 0.25, and TE 
polarization is considered. The Bragg- and Mie-gaps obtained from the band structure calculations are marked by dark green, 
(b, c) The photonic band structure for the 2D square lattice of rods with £ = 60 for TE- and TM-polarization respectively for 
the r —»• X scan of the wave vector k. 



FIG. 3: Experimental setup for transmission mea¬ 
surements and schematic of the lattice transforma¬ 
tion. An anechoic chamber with a metacrystal composed of 
10 tubes in length (x-axis) and 5 tubes in width (y- axis), 50 
tubes in total. Each tube has a length of 1 m (along 2 -axis), 
an outer radius of 2 cm and an inner radius of 1.7 cm. All 
tubes are connected in a closed series circuit together with a 
water heater and a temperature stabilizer. Insert: Schematic 
illustration of the transformation of the 2D square lattice com¬ 
posed of dielectric rods as the lattice constant a increase. The 
scheme shows the square lattice from above. The rods with 
dielectric permittivity € and fixed radius r are embedded in 
air (Eair = 1)- The material is homogeneous along the 2 di¬ 
rection and periodic along x and y. 


vector wave problem can be reduced to two independent 
scalar equations. As a result, the modes of 2D PhCs can 
be classified as either transverse-electric modes TE (E X1 
E y , H z ) with the electric field confined to the x — y plane 
and the magnetic field polarized along the axes of the 
rods, or transverse-magnetic modes TM (H x , H y , E z ) 
for which the magnetic field is oriented perpendicular to 
the rod axis 2 . For the TM polarization such structures 
can not show any magnetic activity^, so we only consider 
TE-polarization. 

Figure 4 demonstrates excellent agreement between 
the measured and calculated transmission spectra both 
for r/a - and e - dependencies. For the r/a-scenario, we 
started our experiment from the PhC phase (r/a = 0.08) 
and compressed the structure by decreasing the lattice 
constant a (Fig. 4a). The lowest Bragg gap demon¬ 
strates remarkable shift towards high-frequency region, 
then nearly almost disappears at r/a « 0.09 when it 
approaches the unchanging TEoi Mie gap from the low- 
frequency side, and finally appears again at r/a ~ 0.095 
on the high-frequency side of the Mie gap. As a result, 
the lowest Bragg band and the TEoi Mie band change 
their positions in the energy scale inaugurating the tran¬ 
sition into MM phase. In contrast, for the ^-scenario, the 
Bragg gap position changes very slowly as opposed to the 
strong shift of the TEoi Mie band that finally leads to 
the same general result (Fig. 4d). 


DISCUSSION 


netic wave incident perpendicular to the tube axis 2 , the 
two transverse polarizations appear decoupled and the 


Figure 5a shows our main result: the PhC-MM phase 
diagram for the 2D square lattice of circular rods. The 
axes represent the dielectric permittivity 6 and filling 
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FIG. 4: Transmission spectra for a square lattice of circular rods: Experiment and numerics. TE polarization, 
a-c r/a-scenario, e = 62. The calculated b and experimentally measured c transmission spectra of the metacrystal as a function 
of the filling ratio r/a. The calculated spectra are shown in the range 0.05^r/a^0.21, the experimental spectra are shown in 
the range 0.08^r/a^0.19. d-f ^-scenario, r/a — 0.08. The calculated e and experimentally measured f transmission spectra 
of the metacrystal as a function of the dielectric permittivity. The calculated spectra are shown in the range 30^5^94, the 
experimental spectra are shown in the range 62^£^94. The spectra are shifted vertically by the constant value. 



FIG. 5: PhC-MM phase diagram, a PhC-MM phase diagram for a square lattice of dielectric circular rods. TE polarization. 
Blue circles indicate the points where the TEoi Mie gap splits from Bragg band becoming the lowest gap in the spectrum (Fig. 
6). The experimentally investigated regions are marked by red. The dielectric permittivities £ of different materials are marked 
by horizontal lines, b The dependence of the order parameter r/ on the dielectric permittivity £ and the filling factor r/a. The 
order parameter is defined as a difference in the energy between the TEoi Mie and Bragg gaps in the MM phase. 


ratio r/a. The phase diagram was obtained based on mittivity region 1^£^100 with the step r/a = 0.01 for 
50 gap maps (such as in Fig. 2) calculated in the per- TE polarization. The diagram is divided into two areas, 
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which represent the PhC and MM phases. The transition 
from PhC-phase to MM-phase can be identified by crit¬ 
ical points which are defined by a special combination 
of permittivity e and filling ratio r/a. The curvature 
of the second dispersion curve around the T point de¬ 
creases and finally becomes perfectly flat along the TX 
path (at 5 = 19.5, r/a = 0.25, Fig. 6d). With future in¬ 
crease in e or change in r/a, the TEoi Mie gap completely 
splits from complicated coupled Mie-Bragg band becom¬ 
ing the lowest gap in the spectrum and the MM-phase 
appears. Within this Mie gap the magnetic permeability 
is negative and all waves in the medium are evanescent. 
Since all the Bragg gaps are located at higher frequen¬ 
cies, any diffraction losses in the vicinity of the Mie reso¬ 
nance are absent. The transition from the PhC phase to 
the MM phase is completed and the medium is classified 
as the MM with /r < 0. A phase transition is usually 
characterized by an order parameter that basically arises 
from symmetry breaking but can also be defined for non¬ 
symmetry-breaking transitions^. Here for the unusual 
phase transition PhC—>MM we define the order parame¬ 
ter r] as a measure of splitting in the energy scale between 
the TEoi Mie and Bragg gaps. It is zero in the PhC phase 
and nonzero in the MM phase. Figure 5b presents the 
order parameter depending on both external variables 5 
and r/a. 

Now we review some features of the PhC-MM phase 
diagram. At low filling ratio r/a < 0.05, the concentra¬ 
tion of the completely isolated rods in the lattice is too 
small for achieving negative average permeability at rea¬ 
sonable values of 5. In contrast, when distance between 
rods become very small (r/a > 0.45), the inter-rod inter¬ 
action become strong and the hopping of the Mie pho¬ 
tons from one rod to another dominates. This means 
that Mie resonances are not localized anymore and they 
finally exhibit the character of propagating Bloch waves. 
All requirements for the appearance of the MM-phase are 
broken and the PhC-phase extends its region of existence 
to very high 6. 

Identification of candidate dielectric low-loss materials 
is critical in designing a practical MM for different ranges 
of electromagnetic spectrum. As an example, we plot in 
Fig. 5a the regions of the MM-phase for water in mi¬ 
crowave range (region 55 ^£^87)22, chalcogenide glasses 
in IR (e « 41)22 and tellurium in visible (e ~ 33). The 
2D square lattice of silicon (e « 16)22 rods has only PhC- 
phase. 


Conclusions 

In this paper we integrate a 2D square structure com¬ 
posed of dielectric circular rods within a general phase di¬ 
agram in the permittivity-filling ratio plane. The bound¬ 
aries of the MM-phase were obtained theoretically and 
confirmed experimentally. Based on the proposed ap¬ 
proach one can obtain different PhC-MM phase diagrams 
by altering the dimensions, symmetry, composition, size 


and geometry of structural elements within a cell. The 
underlying mechanism of the PhC-MM transition pro¬ 
vide unique criterion for goal-oriented search of novel low 
- loss MM with a wide range of applications at infrared 
and visible frequencies. 
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METHODS 

The scattering properties of an isolated dielectric cir¬ 
cular rod and 2D square lattices composed from such 
rods were calculated using three independent methods. 
We calculated (i) spectra of the Mie scattering by an iso¬ 
lated rod; (ii) photonic band structure of a 2D square 
lattice composed of rods; and (iii) transmittance of a 2D 
square lattice composed of 10 rods in length (along x- 
axis). All data sets calculated for a wide range of the rod 
permittivity 1^£^100 with the step of A£ = 1 and in the 
range of filling ratio 0^r/a^0.5 with the step of 0.01. 

First we consider the Mie scattering by an isolated rod. 
The far-held scattering can be expanded into orthogonal 
electromagnetic dipolar and multipolar modes, described 
by circular Lorenz-Mie resonant coefficients a n (TE po¬ 
larization) and b n (TM polarization) corresponding to 
electric and magnetic moments, respectively^. Figure 
6 presents the numerically calculated spectra of the Mie 
scattering efficiency Q SC a,n = f l a n | 2 for the dipole TEo 
and multipole TE n (n^l) modes in the range of di¬ 
mensionless frequencies a/X from 0 to 0.8. Each mode 
consists of an infinite set of quasi-equidistant harmonics 
TE n £ (fc^l). The calculations reveal a strong decrease of 
the Mie eigenfrequencies and narrowing of the resonant 
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FIG. 6: (Left-hand sections in all panels) Calculated Mie scattering efficiency Q S ca,n for an isolated dielectric circular rod for 
TE n (n > 0) modes. (Central sections) The band structure for 2D square lattice of rods with r = 0.25a in air ( e a i r = 1) for the 
TE polarization. The band structure is shown between the T point (wave vector k = 0) and the X point (\k\ = tv/ a along the 
x direction). (Right-hand panels) The transmittance calculated for 10 lattice layers of the 2D square structure of rods in air 
for the TE polarization. The frequency and wave vector are plotted in dimensionless units a/ A, where A is vacuum wavelength 
and a denotes the lattice constant. a £ = 4, b £ = 9, c e — 12.8, d e = 19.5, e e — 25, f £ = 40. 


Mie bands with e increasing, as can be clearly seen for 
the lowest TEoi magnetic Mie resonance in Fig. 6. 

The photonic band structures for the 2D square lat¬ 
tice of circular rods were computed numerically using 
the plane wave expansion method^. The square 2D lat¬ 
tice has a square Brillouin zone, with the three special 
points T, X, and M corresponding to k « 0, k = 
and k « respectively. We restrict here our dis¬ 

cussion to the case for which the wave vector k of the 
eigenmodes is oriented in the T aX direction of the first 
Brillouin zone displaying the lowest Bragg frequencies 


in comparison with the T —>-M direction (Fig. 6). The 
photonic band structures were obtained with 128 by 128 
plane waves. 

The transmission spectra were calculated by using the 
CST Microwave Studio software for the wave vector of the 
incident beam parallel to the T —>>X direction. We con¬ 
sidered the 2D square lattice with different thicknesses 
along x axis which leads to creation of the photonic band 
gaps, including thicknesses of 10 lattice layers as shown 
in Fig. 6. 
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